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Majorana braiding with thermal noise 
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We investigate the self-correcting properties of a network of Majorana wires, in the form of 
a trijunction, in contact with a parity-preserving thermal environment. As opposed to the case 
where Majorana bound states (MBSs) are immobile, braiding MBSs within a trijunction introduces 
dangerous error processes that we identify. Such errors prevent the lifetime of the memory from 
increasing with the size of the system. We confirm our predictions with Monte Carlo simulations. 
Our findings put a restriction on the degree of self-correction of this specific quantum computing 
architecture. 
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Introduction - Taking advantage of topological states of 
matter to encode and process quantum information is 
viewed as a promising route towards quantum compu¬ 
tation. The idea is that the quantum information is 
naturally protected from decoherence so that quantum 
gates can be performed reliably by braiding anyonic ex¬ 
citations. Recent work has focused on the realization of 
the topological phase as described by the Kitaev wire 
model mm- Unpaired Majorana modes are predicted 
in this model, and, under braiding in a network of ID 
wires, these modes will behave as Ising anyons, enabling 
a form of topological quantum computation. Intensive 
theoretical nu and experimental mm investigations 
of nanowire hybrid systems have been undertaken to es¬ 
tablish the existence of these “Majorana particles”, with 
the hope that a network of branched nanowires pro¬ 
vides a suitable geometry for the moving of these par¬ 
ticles 113 , 12 ]. 

While topology does give some measure of robustness, 
especially in the face of static disorder H 33 QH, it has been 
noted that Majorana qubits are not entirely immune to 
decoherence. For example, these qubits are dephased by 
the induced splitting of nearby Majorana states due to 
tunneling m- It is also clear that it is essential to avoid 
parity-changing excitations, that is, injection of individ¬ 
ual quasiparticles from the environment, which are im¬ 
mediately destructive of Major ana-qubit coherence US¬ 
ED! . The decoherence of Majorana qubits induced by 
thermal excitations above the superconducting gap has 
been investigated in Refs. I 2 T 1 and l22l Additional sources 
of noise in a trijunction setup [ 23 ] and when Majorana 
bound states (MBSs) are physically moved [ 24 H 281 have 
also been recently studied. 

Starting from a microscopic model, here we perform a 
detailed analysis of the effect of a parity-preserving en¬ 
vironment on the fidelity of the topological qubit. Our 
main focus is the interplay between thermal errors and 
the adiabatic motion of MBSs, within the trijunction 
scheme of M- Our calculations here, using a standard 
treatment of a generic local bosonic environment at low 
temperature, show that the parity-preserving environ¬ 
ment produces unavoidable dangerous errors arising from 


the motion of MBSs. These errors put an upper limit on 
the lifetime of the qubits, which cannot be improved by 
keeping the MBSs further apart. The motion of MBSs 
causes a spatial separation of thermally produced quasi¬ 
particle pairs, such that they act to effectively change the 
parity, just as parity-breaking environments do. 

Kitaev model- The Kitaev wire [2], an archetypical exam¬ 
ple of a system supporting Ising anyons, consists of a ID 
chain of hopping fermions coupled by a superconducting 
pairing term, 

L L -1 

Hs(t) = - ^/ij( T )aja r ^](ta]ffl i+ i-Aa j a i+ i+h.c.). 

3=1 3 =1 

(i) 

The first term describes a site- and time-dependent chem¬ 
ical potential Rj(r) < 0 . The second and third terms de¬ 
scribe respectively nearest-neighbor hopping with t > 0 
and superconducting pairing with A = |A|e^. 

We first analyze on the case fij = 0 \/j and express 
Hs in terms of Majorana-mode operators aj = (72,7-1 + 
272.7 )/2 which satisfy {77,7/c} = 2 Sjk- For the parameters 
we study in this paper, t = |A| and 6 = 0 , Hs reduces to 

l- 1 

H s = -A X H2J+172 j ■ (2) 

i =1 

Note that the end Majorana modes, 71 and 72 l, are de¬ 
coupled from the rest of the chain; thus, a zero-energy 
delocalized fermionic mode exists with creation operator 
= (71 + ^72 l)/2. Hs can be fully diagonalized using 
the eigenoperators dj = (727 + 2727+1 )/2 and written as 

"> = AEb,'( 24 /,- 1). 

For jij 7^ 0 , MBSs localized near the end of the wire 
persist as long as \fij\ < 2£; this is called the topological 
phase. For \fij\ > 2 1 the localized modes disappear and 
we enter the nontopological phase. Deep in the nontopo- 
logical phase, with \fij\ t, A, the Majorana Hamilto¬ 
nian approaches H^ ontop = —ifi/2 727-1727 and the 
Majorana modes are paired in a shifted way as compared 
to Eq. see Fig. Ek 

More generally, MBSs occur at the junction between 
topological and nontopological segments of the wire. By 
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FIG. 1. a) Representation of the Kitaev wire supporting a 
nontopological (gray) and a topological (black) segment. The 
large black dots represent the fermionic sites of Eq. |l]) . The 
smaller dots below represent the Majorana modes of Eq. §. 
The lines between the small dots represent the coupling be¬ 
tween Majorana modes in the limit \/ij\ t, A (for the non¬ 
topological segment) and for the values fij = 0 and t — A (for 
the topological segment). The excitations above the ground 
states are local and can be understood as ip particles as de¬ 
picted. b) Space profile of the chemical potential pj corre¬ 
sponding to the situation in a). The chemical potential in 
the nontopological segment is made non-uniform in order to 
localize the T/Z-paxtides. c) Configuration of ip particles cor¬ 
responding to logical A, Y, and Z Pauli errors in the trijunc¬ 
tion. 


varying the chemical potential, one can increase or de¬ 
crease the size of the nontopological segments and thus 
move the position of the MBSs. As proposed in Ref. na, 
this technique can be used to braid MBSs in a nanowire 
trijunction setup (Fig. [ 2 | and perform gates in a topo¬ 
logically protected fashion. In order to stay within the 
ground-state subspace, the motion of the MBSs must be 
adiabatic. This means that the chemical potentials Pj(r) 
must be varied on a timescale slow compared to 1/A. 
Specifically, here we vary the chemical potentials linearly 
with time. 

Encoding / Error Correction - For the sake of clarity, we 
first focus on the encoding of a logical qubit in a single 
Kitaev wire with pj = 0 and then we generalize the dis¬ 
cussion to the branched wire forming a trijunction. The 
ground-state subspace of Hs in Eq. © forms a stabilizer 
code with stabilizer operators Sj = [HE 29 ] : 

that is, the logical qubit states |0) and |1) are invariant 
under the action of the Sj operators with j = 1 , ..., L— 1 . 
Furthermore, they are defined by the phase they acquire 
under the action of the parity operator So = 271721,, 
namely So\0) = | 0 ) and So\l) = Hi). 

In this context, excitations above the ground states 
are interpreted as particles; whenever a stabilizer is vi¬ 
olated, then we say that the wire supports a ip parti¬ 
cle at the position of the violated stabilizer, see Fig. Eh 
Also here we represent a parity flip So —So by draw- 



FIG. 2. a) Trijunction of size L — 7, composed of horizon¬ 
tal and vertical wires coupled by hopping and superconduct¬ 
ing pairing. The zero-energy MBSs 71 , 2 , 3,4 are shown, b) 
Braiding sequence considered in this work; MBSs 1 and 3 
are exchanged, while MBSs 2 and 4 remain immobile, c), 
d ): Sequence of two equiprobable error processes. The two 
sequences lead to the exact same error syndrome. After ap¬ 
plication of our error correcting algorithm to the evolution at 
v), the outcome in c) is faulty while the one in d) is successful. 


ing a ip particle inside one of the MBS. In the ground- 
state subspace, drawing the ip inside the MBS 71 is in¬ 
terpreted as a logical X flip of the qubit state, drawing 
it inside the MBS 72 l corresponds to a logical Y flip, 
while drawing it inside both of them corresponds to a 
logical phase flip XY oc Z = So- Note that here we con¬ 
sider solely parity-conserving perturbations and ip parti¬ 
cles are thus created in pairs. A string operator creating 
pairs of ip’s corresponding to, say, Sj = —1 and Sk = — 1 
is Sjk = 72j+i72j+2 • • • 72fc. String operators give us a rig¬ 
orous way to define fusion of the different excitations, and 
thus perform error correction; two ip excitations S a = — 1 
and Sb = — 1 are fused and annihilated by applying S a b- 
If one wants to fuse a ip excitation S a = — 1 with the MBS 
71, one should apply So a 5 this moves the ip from the bulk 
to inside the MBS, thereby flipping the parity So — So • 
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This phenomenology corresponds to the standard fusion 
rules for Ising anyons m with topological charges cr, ip, 
and 1 : ax<j = l + ip,ipxip = l, and a x ip = a. MBSs 
are identified with the Ising anyons a, because they fol¬ 
low the prescribed fusion rules (as just shown) as well as 
braiding rules m- 

We now discuss the trijunction setup. We consider a 
scenario with four MBSs 71-4, using a scheme in which 
71 and 73 are braided while 72 and 74 remain immobile, 
see Fig. ©>. In such a case, the ground-state subspace is 
4 -fold degenerate and we follow the procedure presented 
in Ref. [31 to encode the logical qubit in a fixed-parity 
sector, say 2717227374 = + 1 . Similar to above, we have 
2717210) = 27374|0) = |0) and 1717211) = 27374II) = -|1>. 
Again, excitations are localized and can be understood 
in terms of ip particles. Logical X-, Y-, and Z-errors are 
present when ip particles reside inside MBSs, see Fig. [TJ^. 

Pairs of ip particles can be created, annihilated or can 
hop due to interaction with a bath. But we presume the 
possibility of final error correction that undoes the effect 
of the thermal environment and protects the stored quan¬ 
tum information. We thus assume that it is possible to 
identify the position of the ip particles in the bulk of the 
trijunction setup and thus identify the error syndrome. 
Note that the fusion rule a x ip = a indicates that a ip 
inside the MBS is invisible to error correction. Having 
the syndrome information in hand, one performs a set 
of unitary operations having the effect of annihilating ip 
particles in order to map back into the ground-state sub¬ 
space and hopefully recover the stored qubit. The simple 
algorithm we use is an adaptation of the algorithm pre¬ 
sented in Ref. [221 US- The main idea is to fuse pairs of 
nearby particles (MBS or ip). 

Coupling to a Thermal Bath- We see that, over time, 
the fundamental ip excitations cause the destruction of 
the encoded information. Here we derive a master equa¬ 
tion for the trijunction when such excitations result from 
coupling to a thermal bath. We start with the following 
microscopic model: 

H(t) = Hl rii (r) + H b + H sb • ( 3 ) 

Here is the adiabatically varied trijunction 

Hamiltonian, consisting of hopping and pairing terms 
as in Eq. ( 1 ) for each segment of the trijunction, Hb 
is a bosonic bath Hamiltonian and Hsb = ® 

( 2 OjOj — 1 ) = —i JA Bj ® 72 j—i 72 j is a conventional 
system-bath coupling. Importantly, Hsb conserves the 
overall parity of the system; system operators 72j-i72j 
anticommute with two stabilizers Sj-i and Sj, thus ip 
excitations are created in pairs. 

Following the Davies prescription [ 33 ], we obtain the 
master equation, in Lindblad form, for the system density 
matrix in the adiabatic limit, 

p s (r) = -i [-Hs 1 J (r),p s (r)] +V{p s (r )), ( 4 ) 


with 

^(Ps(r)) = E (^(w(r))p5(r)A J (w(r)) t ) 

i,3 uj(t) 

-\{ A3 ( w ( r )) fA ‘( w W). AsO)}^ , (5) 

where lj(t) describes the different energy gaps at time 
r between the different eigenbranches of iL^ riJ (r). The 
Davies approach gives jump operators A 1 (uj(t)) satisfy¬ 
ing detailed balance, that cause transitions between dif¬ 
ferent system eigenbranches separated by energy c o(r) at 
time r. A crucial property of our system Hamiltonian is 
that the jump operators remain strictly local [ 33 ] . 

The diagonal matrix elements of ps(t) decouple from 
the off-diagonal elements [ 33 ] so that one obtains a Pauli 
master equation for the population in an instantaneous 
eigenstate |n(r)) of iL^(r), 

——VF(n(r)|m(r))P(m(r), r) — m o n , 

m(r) 

(6) 

W(n(r)\m(r)) = 7(w mn (r))|(TO(r)|H* m ’ l (w mn (r))|n(r))| 2 

and 7 (cj) = k |cj/( 1 — exp(—/ 3 a;))| is the bath spectral 
function with inverse temperature /? and with coupling 
constant k. Note that 7(0;) does not depend on posi¬ 
tion since we assume the baths coupled to each site to 
have identical form and thus to have the same spectral 
function. 

For the sake of illustration, we first consider here the 
time-independent case, i.e. when the MBSs are immobile, 
and identify which transitions are caused by the jump 
operators. In the case of moving MBSs many more fun¬ 
damental processes are allowed, but we believe that the 
following list gives some intuition for the error processes 
that also occur in the more general case, i) Creation 
or annihilation of a i^-pair in the bulk, with energy cost 
± 4 A. ii) Hopping of a ip to a nearest neighbor site in the 
bulk, with energy cost 0. iii) ^-pair creation (annihila¬ 
tion) at a boundary supporting a MBS, with energy cost 
± 2 A. iv) Hopping of a ip into or out of a neighboring 
MBS, with energy cost ±2A. In Ref. [ 33 ], we present an 
exhaustive list of the hundreds of allowed error processes 
that occur in the trijunction during adiabatic braiding. 

It is worth pointing out that the system-bath interac¬ 
tion in Eq. ©> does not support the creation of ip par¬ 
ticles in the nontopological segments of the trijunction. 
Indeed, as the chemical potential becomes very negative, 
the eigenstates become very close to those of iL^ ontop and 
[iL^ ontop , Hsb] = 0. However, this does not mean that no 
ip particles will ever be present in the nontopological seg¬ 
ments. To understand that, one needs to investigate the 
interaction between moving MBSs and ip particles. We 
note: i) When a MBS moves over an existing ip particle, 
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the ip particle is transferred from the topological seg¬ 
ment into the nontopological one. For clarity we use the 
notation ip' for an excitation in the nontopological seg¬ 
ments. In the reverse scenario, when a MBS moves over 
a ip ', then the excitation transfers back to the topological 
segment, becoming a ip. ii) A ip particle inside a MBS 
remains inside the MBS during the motion. These two 
properties can be proved straightforwardly by diagonal- 
ization of a four-site model [33]. In our error-correcting 
algorithm presented above, ip and ip' particles are treated 
in the exactly same way. In particular, during error cor¬ 
rection it is assumed that the position of excitations can 
be measured everywhere, not only in the topological seg¬ 
ments. 

Dangerous Error Processes- We show that some elemen¬ 
tary error processes lead to a lifetime of the stored logical 
qubit that cannot be improved by increasing the size of 
the trijunction setup. This is in contrast with the case of 
immobile MBSs, where the lifetime grows with the sys¬ 
tem size [0 3 33] ; in Ref. [33] we show that the lifetime 
grows logarithmically with L in the high-temperature 
limit. Immobile MBSs thus represent a good quantum 
memory. However, a logarithmic scaling of the lifetime 
is generally too weak if braiding is performed since the 
time to adiabatically braid the MBSs scales linearly with 
L [34]. Here we show that the lifetime does not actu¬ 
ally scale at all with L when MBSs are moved because 
of dangerous error processes. These error processes thus 
put a restriction on the degree of self-correction of this 
specific quantum computing architecture when braiding 
is executed. 

Consider the creation of a V’-pair, one inside a MBS 
and one in the bulk of the trijunction. Since a ip inside a 
MBS cannot escape when the MBS is moved, this ip par¬ 
ticle will be dragged along during the braiding motion. 
In other words, an originally local error process becomes 
highly nonlocal due to the braiding motion. For example, 
the error sequence depicted in Fig. leads to a logical 
X error after our error correcting algorithm has been ap¬ 
plied, although only two fundamental error processes oc¬ 
curred. Dangerous error processes will generally lead to 
X and Y logical errors and our error correcting algorithm 
will fail independent of the system size. It may appear 
surprising that X and Y errors are at all possible since 
the total parity of the system is conserved. However, this 
can happen because the parity of individual topological 
segments is not conserved since MBSs are moved and ip 
particles can be transferred to nontopological segments 
as mentioned above. 

One natural question is whether a better algorithm 
would be able to take into account the non-locality in¬ 
troduced by braiding. Unfortunately, this is impossible, 
as shown by Figs. and d. Here we see two distinct 
sequences of two fundamental error processes leading to 
exactly the same error syndrome. The difference between 
the two final states is the occurrence of the ip particles 



FIG. 3. a) Probabilities px+Y (solid) and pz (dashed) as 
a function of time for trijunctions with L — 19 , 29 , 49 . The 
elementary time step to move a MBS to a nearest-neighbor 
site is 10 4 /A. The braiding time is Tbraid = (2 L + 2)10 4 /A; 
a dashed vertical line specifies when Tbraid is achieved. For 
times r < Tbraid, we execute braiding motion incompletely 
while the coupling to the external bath is on. We then fin¬ 
ish the braiding motion unitarily without any coupling to the 
thermal bath, b ) Artistic representation of the total proba¬ 
bility of failure as a function of time for a trijunction setup 
of a given length where three braiding motions, separated by 
a time interval Sr, are executed. If the coupling to the ex¬ 
ternal bath is small, the only contribution to ^failure is due to 
dangerous errors during braiding. 

inside MBSs. But these ip ’s are invisible to error cor¬ 
rection, so these two situations produce the same error 
syndrome and are not distinguishable by any error cor¬ 
recting algorithm. If one of these cases is successfully 
recovered by error correction, the other will lead to fail¬ 
ure. Our discussion is independent of the size of the 
trijunction; this implies that the lifetime of the stored 
information will not grow with system size. Our study 
thus provides evidence that active error correction dur¬ 
ing braiding would be necessary in topological quantum 
computation [351438] . 

Monte Carlo- We confirm our predictions by simulating 
the Pauli master equation § with standard Monte Carlo 
methods. We have calculated the rates of all the allowed 
error processes by numerically diagonalizing six-site tri¬ 
junctions, see Ref. [35] for more details. We focus on the 
low-temperature regime /3 = 4/A. A discussion of the 
high-temperature case can be found in Ref. [33] . 

Fig. [ 3)2 shows the probabilities px+Y , and pz of X and 
y, as well as Z logical errors as a function of time for tri¬ 
junctions with different L [39 . A single braid (Fig. 2b) 
is executed at the beginning and the MBSs are then left 
immobile for the remaining time. px+Y increases signif¬ 
icantly during the braiding period. This is in agreement 
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with the above discussion of dangerous processes that oc¬ 
cur solely during the motion of MBSs. Also, the growth 
of px+Y at small times is found to be independent of 
L, confirming that dangerous processes lead to a lifetime 
independent of the system size. On the other hand, Z 
logical errors are due either to a high density of parti¬ 
cles or to diffusion of a ^ particle over a distance greater 
than L/ 2. In both cases error correction will eventually 
fail due to an accumulation of error processes. There¬ 
fore, pz stays small during the (short) braiding period 
and increases at larger times. 

Conclusion- While our results put restrictions on the 
quantum coherence that may be expected in one par¬ 
ticular braiding scheme, there are many other scenar¬ 
ios in which MBSs could have unique strengths in the 
processing of quantum information. For example, the 
present analysis does not restrict the feasibility of the 
interaction-induced braiding concept of Refs. Haim ; we 
hope to pursue an extension of our system-environment 
treatment to this case. 
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